A new numerical method based on LZB is proposed for nonsmooth dynamics on astronautics separation device in guided stage. Using Euler-Lagrange equations and the LZB method, the general calculating framework describing the whole dynamical motion including contact, collision with friction is established. Compared with FEM simulation time, numerical simulation implies that this method is valid.